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We give a statistical characterization of states with nonzero 
winding number in the Phase Turbulence (PT) regime of the 
one-dimensional Complex Ginzburg-Landau equation. We 
find that states with winding number larger than a critical 
one are unstable, in the sense that they decay to states with 
smaller winding number. The transition from Phase to De- 
fect Turbulence is interpreted as an ergodicity breaking tran- 
sition which occurs when the range of stable winding num- 
bers vanishes. Asymptotically stable states which are not 
spatio-temporally chaotic are described within the PT regime 
of nonzero winding number. 

PACS: 05.45.+b,82.40.Bj,05.70.Ln 



Spatio-temporal complex dynamics jl],^| is one of the 
present focus of research in nonlinear phenomena. Much 
effort has been devoted to the characterization of differ- 
ent dynamical phases and transitions between them for 
model equations such as the Complex Ginzburg-Landau 
Equation (CGLE) One of the main questions 

driving these studies is whether concepts brought from 
statistical mechanics can be useful for describing com- 
plex nonequilibrium systems . In this paper we give 
a characterization of the spatio-temporal configurations 
that occur in the Phase Turbulence (PT) regime of the 
CGLE (described below), for a finite system, in terms 
of a global wavenumber. This quantity plays the role of 
an order parameter classifying different phases. We show 
that in the PT regime there is an instability such that a 
conservation law for the global wavenumber occurs only 
for wavenumbers within a finite range that depends on 
the point in parameter space. Our study is statistical 
in the sense that averages over ensembles of initial con- 
ditions are used. Our results allow a characterization 
of the transition from PT to Defect or Amplitude Tur- 
bulence (DT) (another known dynamical regime of the 
CGLE) as the line in parameter space in which the range 
of conserved global wavenumbers shrinks to zero. 

The CGLE is an amplitude equation for a complex field 
A(x, t) describing universal features of the dynamics of 
extended systems near a Hopf bifurcation 



[Q, chemical turbulence p5[ , bluff body wakes Jl6| , 
among other systems can be described by the CGLE in 
the appropriate parameter range. We will restrict our- 
selves in this paper to the one-dimensional case, that is 
A = A(x,t), with x € [0, L]. For this situation a ma- 
jor step towards the analysis of phases and phase transi- 
tions in (|]) was the identification (3[jj| of different chaotic 
regimes in different regions of the parameter space [c\ , C2] 
(see Figure 1). 
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FIG. 1. Regions of the parameter [ci , C2]-space for the 
CGLE displaying different kinds of regular and chaotic be- 
havior. Lines Li,Ls were determined in 
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Binary fluid convection |13] , transversally extended lasers 



Eq. ([!]) has plane-wave solutions Ak = \/l — k 2 e tkx 
with k e [—1,1]. When c\C2 > —1 there is a range of 
wavenumbers [— fc^fe^] such that the plane- wave solu- 
tions with wavenumber in this range are linearly stable. 
They become unstable outside this range (the Eckhaus 
instability Q). The limit of this range, ks approaches 
zero as the product c\C2 approaches —1, so that the range 
of stable plane waves vanishes by approaching from be- 
low the line C1C2 = — 1 (the Benjamin- Feir or Newell line, 
labeled BF in Fig.l). Above that line no plane wave is 
stable and different turbulent states exist. The authors 
of (^-^| identified three different regimes in different re- 
gions above the BF line (Fig.l): PT, DT, and bichaos. 
Among these regimes, the transition between PT and 
DT has received special attention P JlQ| , p7| . In spite of 
the fact that there are some indications that this transi- 
tion can be ill-defined in the L — * 00 limit |@,0,Eo[, the 
PT regime is robustly observed for any finite size system 
and for finite observation times, with the transition to 
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DT appearing at a quite well defined line (Li in Fig.l) 
@. In the DT region the modulus |A| of A = \A\e^ be- 
comes zero at some instants and places (called defects), 
so that the phase (f> becomes undefined and the winding 
number v = ^_ Q x (j)dx changes value during evolu- 
tion. In contrast, dynamics maintains the modulus of A 
far from zero in the PT region, so that v is thought to be 
a constant of motion there. A global wave number of the 
configuration can be defined as k = 2irv/L. These differ- 
ent regimes were originally identified from the analysis of 
the space-time density of defects. If this picture is cor- 
rect, one can speculate that the transition between DT 
and PT would be a kind of ergodicity breaking transition 
p8[ as in other systems described by statistical mechan- 
ics. DT would correspond to a "disordered" phase and 
v classifies different "ordered" phases in PT. However, 
we note that most studies of the PT regime have only 
considered in detail the case of v = 0. In fact the phase 
diagram in Fig.l was constructed for this case. In order 
to provide a better understanding of the PT-DT transi- 
tion we undertake in this Letter a systematic study of 
PT configurations with v 7^ . 

Typical configurations of the PT state of zero winding 
number consist of pulses in \A\, corresponding to phase 
sinks, that travel and collide rather irregularly on top of a 
k = unstable background wave (that is, a uniform oscil- 
lation) HH . The phase of these configurations strongly 
resembles solutions of the Kuramoto-Shivashinsky (KS) 
equation. Quantitative agreement has been found be- 
tween the v = PT states of the CGLE and solutions of 
the KS equation near the BF line j[o|. The more obvi- 
ous effect of a non-zero v is the appearance of a uniform 
drift added to the irregular motion of the pulses. In ad- 
dition Chate reported an earlier breakdown of the 
PT regime when v 7^ 0. Our results below show that not 
all the winding numbers are in fact allowed in the PT re- 
gion at long times. PT states with too large \v\ are only 
transients and decay to states within a band of allowed 
winding numbers. The width of this band shrinks to zero 
when approaching the line L\. In addition we find that 
the allowed non-zero winding number states are not of 
a single type. We have identified three basic types of 
asymptotic states for v 7^ 0. Two of them are different 
from the usually described PT states in the sense that 
they do not exhibit spatio-temporally chaos. 

In order to study the dynamics of states with v 7^ 
we have performed simulations extensively covering the 
PT region of parameters of Fig. 1. Only a small part 
of the simulations is shown here, and the rest will be 
reported elsewhere. We use a pseudo-spectral code with 
periodic boundary conditions and second-order accuracy 
in time. Spatial resolution was typically 512 modes, with 
runs of up to 4096 modes to confirm the results. We 
work at fixed system size L = 512. The initial condition 
in our simulation is a plane wave of the desired winding 



number, slightly perturbed by a white Gaussian random 
field. The initial evolution of the spatial power spectrum 
is well described by the linear stability analysis around 
the initial plane wave: typically the perturbation grows 
mostly around the most unstable wavenumbers identified 
from such linear analysis. After some time the system 
reaches a state similar to the v = PT, except for a non- 
zero average velocity of the chaotically travelling pulses. 
We call this state riding PT. Its spatial power spectrum is 
broad and unsteady, with the more active wavenumbers 
located around the one determined by the initial winding 
number. We observe that when this winding number 
is small, it remains constant in time, and the system 
either remains in the riding PT state or approaches one of 
the more regular asymptotic states that will be described 
below. If \v\ is initially too high, the competition between 
wavenumbers leads to phase slips that reduce \i>\ until a 
value inside an allowed range is reached. Then the system 
evolves as before. 

We present in Fig. 2 the temporal evolution of P(t), 
the average of v{t) over 50 independent realizations of 
the random perturbation added to the initial plane wave 
for a fixed point in parameter space. The variance among 
the sample of 50 realizations is also shown. Three initial 
values Vi of the winding number are shown. v(t) typically 
presents a decay from i/i to the final winding number Vf. 
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FIG. 2. a) Temporal evolution of v(i) for three different ini- 
tial winding numbers Vi = 19 (solid), 15(dotted), 10 (dashed) 
. ci = 2.1, C2 = —0.75. b) Winding number standard devia- 
tion a. 

The decay is found to take place in a characteristic 
time t that we quantify as the time for which half of the 
jump in v has been attained. Figure 3 shows 1/r for 
different values of i/j. The different curves correspond to 
different values of C2 with fixed c\. Similar results were 
obtained for C2 fixed and varying c\. r increases with an 
apparent divergence as vi approaches a particular value 
v c which is a function of C\ and C2. We estimate this v c 
by fitting linearly the data for 1/r. Other fits involving 
non-trivial critical exponents have been tried, but they 
do not improve the simpler linear one in a significant 
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manner. A very similar value of v c is obtained by simply 
determining the value of Vi below which v(t) does not 
change in any of the realizations. Values of v c from some 
of the simulations are in the inset of Fig. 3. v c vanishes 
as C2 approaches the transition line L\ (or L% when pass- 
ing through the bichaos region). For example, the linear 
fitting of the data in the inset of Fig. 3 and extrapola- 
tion towards zero v c reproduces the value for L\ of |}]|| 
(c2 ~ —0.9 for ci = 2.1) within the fitting error in C2 of 
±0.02. 
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FIG. 3. Inverse of the characteristic time for winding num- 
ber relaxation as a function of the initial winding number. 
The value of ci is fixed (ci = 2.1) and c 2 varies from near 
the B-F line (c 2 = -1/2.1) to the Li line (c 2 ~ -0.9). Dif- 
ferent symbols correspond to c 2 = —0.6 (+), c 2 = —0.7 (*), 
c 2 = -0.75 (O), c 2 = -0.8 (A), c 2 = -0.83 (□). The inset 
shows the critical winding number (y c ) as a function of c 2 . 

The winding number instability found here in the PT 
region is strikingly similar to the Eckhaus instability of 
travelling waves below the BF line of Fig. 1 [§: There 
is a range of allowed winding numbers such that con- 
figurations outside this range undergo phase slips until 
an allowed v is reached. The difference is that below 
the BF line, the attractor for each stable v is a trav- 
elling plane wave of wavenumber k, whereas each i>, or 
equivalent global wave number, characterizes phase tur- 
bulent attractors above the BF line. The allowed range of 
travelling waves shrinks to zero when (ci, ci) approaches 
the BF from below, whereas above BF, the allowed v 
range shrinks to zero when approaching the Li line from 
the right. In this picture, the transition PT-DT appears 
as the BF line associated to an Eckhaus-like instability 
for phase turbulent waves. Since the presence of chaotic 
fluctuations is a kind of self generated "noise" present 
in the system, this winding number instability should 
be compared to the Eckhaus instability in the presence 
of stochastic noise fill, rather than to the determinis- 
tic one. The comparison shows qualitative similarities 
between both cases, but quantitative agreement such as 
similar critical exponents or scaling laws has not been 
found. The comparison is also instructive because it can 



be shown that, for the one-dimensional stochastic case, 
there is no true long range order, and therefore no true 
phase transition in the infinite size limit p0| |. But for 
finite sizes and finite observation times well defined ef- 
fective transitions and even critical exponents can be in- 
troduced @. The PT-DT transition in the CGLE can 
be an effective transition of this kind. In order to further 
characterize the robustness of the effective transition an 
analysis of system size effects should be performed. Pre- 
liminary results indicate that the v c obtained for each 
{c\,Oi) point grows linearly with system size L, as it 
should happen for a well defined extensive quantity. 

Finally we consider the nature of the asymptotic states 
allowed within the band of "stable" v. 
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FIG. 4. Spatiotemporal evolution of d x <j>(x,t) with time 
running upwards and x in the horizontal direction. The 
lighter grey correspond to the maximum value of d x <j>{x, t) 
and darker to the minimum. Different scales of grey are used 
in each case in order to see the significant structures, a) Last 
10 2 time units of a run 10 4 time units long for a tiding PT 
state at ci = 2.1 , and c 2 = —0.83. The initial condition was 
a TW with Vi = 20 that decayed to Vf — —1 after a short 
time, b) Last 10 2 time units of a run 10 5 time units long 
for a quasiperiodic state. The initial condition is random 
noise with an amplitude of 0.05. a = 2.0 , and c 2 = —0.8. 
c)Last 10 2 time units of a run 10 4 time units long for a frozen 
turbulence state. The initial condition is a TW of Vi = 12 
that decayed to Vj = 6 after a short time, a = 1.75 , and 
c 2 = —0.8. d) Spatial power spectrum S(q) as a function of 
wavenumber for the frozen turbulence configuration shown in 
c). This spectrum is constant in time. 
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We have numerically found three basic types of states 
in the PT region of parameters with non-zero v. Fig. 
4 shows in gray levels the value of d x (f)(x, t) as a func- 
tion of x and t. The state shown in the top left is the 
familiar Q riding PT, which is similar to the PT usu- 
ally seen for v = (wiggling pulses in the gradient of 
the phase) except for a systematic drift in a direction 
determined by v. The other two states do not show 
spatio-temporal chaos. They can be described as the 
motion in time of a spatially rigid pattern on the top of 
a plane wave (with k ^ 0) background and with peri- 
odic boundary conditions. The state shown in the top 
right consists of equidistant pulses travelling uniformly. 
They are the quasiperiodic states described in ](|. The 
state shown in the bottom left, that we call frozen tur- 
bulence consists of pulses uniformly travelling on a plane 
wave background, as in the quasiperiodic case, but now 
the pulses are not equidistant from each other. The spa- 
tial power spectrum is shown for this later case. It is a 
broad spectrum in the sense that the inverse of its width, 
which gives a measure of the correlation length, is small 
compared with the system size. This is due to the ir- 
regular positions of the pulses. In addition the spectrum 
is constant in time, which makes this frozen state dif- 
ferent from riding PT. The existence of the two states 
with no spatio-temporal chaos (quasiperiodic and frozen 
turbulence) described above can be understood by ana- 
lyzing the phase equation valid near the BF instability. 
In the case of a non-zero v it contains terms breaking 
the left-right symmetry and it is known as Kawa- 

hara equation |22| . Its uniformly travelling solutions are 
related to the rigidly propagating patterns of Fig. lb) 
and c). These solutions can be analyzed with the tools 
of Shilnikov theory ]23| ]. The details will be discussed 
elsewhere. 

In addition to the pure three basic states there are 
configurations in which they coexist at different places of 
space, giving rise to a kind of intermittent configurations, 
some of them already observed in [[|. The main results 
reported here, that is the existence of an Eckhaus-likc 
instability for phase-turbulent waves, the identification of 
the transition PT-DT with the vanishing of the range of 
stable winding numbers, and the coexistence of different 
kinds of PT attractors should in principle be observed 
in systems for which PT and DT regimes above a Hopf 
bifurcation are known to exist pq| . We note in addition 
that the experimental observation of what seems to be 
an Eckhaus instability for non-regular waves has been 
already reported in p4| . 
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